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, In this work we will study the polarizations of both leptons (r) in the decay channel B — > K*t^t'^ . 

In the case of the dileptonic inclusive decay B K*l , where apart from the polarization 
asymmetries of a single lepton I, one can also observe the polarization asymmetries of both leptons 
simultaneously. If this sort of measurement is possible then we can have, apart from decay rate, 
forward backward asymmetry and the six single lepton polarization asymmetries (three each for l~ 
and t^), nine more double polarization asymmetries. This will give us a very useful tool in more 
' strict testing of the SM and the physics beyond. We discuss the double polarization asymmetries 

^ , of the r leptons in the decay mode B if*r~r"*" within the standard model (SM) and the minimal 

' supersymmetric extension of it. 

o 
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PACS numbers: 13.20He,12.60,-i,13.88+e 



I. INTRODUCTION 



Flavor changing neutral currents (FCNC) in weak decays provide a fertile ground for testing the structure of 
. weak interactions. Since these decays are forbidden in the tree approximation, they go through higher order loop 
effects. Consequently they are sensitive to finer details of the basic interactions responsible for the process and as 
such provide a natural testing ground for any theories beyond the standard model as an example. In the context of 
B-decays, processes involving a dileptonic pair in the final state through the basic quark process b s-^~f ^provides a 
wealth of possible experimental data, accessible in the near future, that can be confronted with theoretical predictions. 
Processes involving this basic quark transition fall into two broad categories, namely the inclusive ones and specific 
exclusive processes. In both these there have been theoretical investigations involving total cross-sections, differential 
cross sections and polarization studies. The last of these, namely polarization studies of the final state particles is 
a particularly useful parameter, since the most popular extension of the standard model (SM) predicts considerable 
modification of their values from SM results EE>I^' Polarizations involving a single lepton have been studied 



extensively in B Xs£~i+R^, B K*£-i+^m,M, B ^ K£~i+ (3, B ^ (7r,p)^-£+ B, ^ i+£-j 
but recently Bensalam et al. [3| have pointed out that the study of simultaneous polarizations of the leptons in the 
final state provides another observable that can be experimentally measured and provides yet another parameter in 
testing models involving physics beyond the standard model. They have, in their, work carried out detailed analysis 
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of the exclusive process B — > Xs£~£'^. On similar double polarization asymmetries of both the leptons this process 
{B Xs£~£~^) should also get major corrections if we consider extension of SM [lo| . 

In Ref. ^ they have confined themselves to the standard model. But as has been emphasized in many works 
[ill ITS that the supersymmetric extension of the SM gives major corrections to the processes based on the quark 
level transitions b — + s£^£'^. Supersymmetry (SUSY) extends the SM list of terms in the effective Hamiltonian and 
associated Wilson coefficients; for the quark level process b s£~£~^it predicts the presence of two new quark bilinears 
in the effective Hamiltonian, namely a scalar and a pseudo-scalar one. These new Wilsons come because of the extra 
Neutral Higgs bosons (NHBs) spectrum of SUSY (and two Higgs doublet model) theories 0, 0, . The effects of 
these new Wilsons on various kinematical variables like branching ratios, lepton pair forward backward asymmetries 
and lepton polarization asymmetries in various inclusive {B Xs£~£+, B Xrf|-£+)JiJlJEi and exclusive 
{B K£-£+,B K*£-£+B £-£+j, B Tr£-£+, B ^ p£-£+etc.) fl B H El 13 13 semi-leptonic and 
pure leptonic {B £-^£+) ;il| decays of B mesons have been studied in great detail. The new Wilson coefficients 
(Cqj and ) are proportional to mimijtan^ l3 and hence can be substantial when the lepton is t and tan (3 is 
sufficiently high. We would like to include the effect of NHBs but at the same time focus on an exclusive process 
B K*£~£~^. Experimentally exclusive processes are easier to study but theoretically involve more uncertainties. 
However for processes like B K*£^ £^the theoretical uncertainties are somewhat in control since the unknown 
hadronic matrix element involved can be related to charged current decay mode of the B meson. The analysis of these 
has been subject to a lot of theoretical attention and one can use the results there as input to theoretical estimates 
for the FCNC process. In this paper we take up the study of this exclusive process for determination of all the three 
polarization parameters, longitudinal, transverse and normal for both the leptons simultaneously. This exclusive 
process is amongst the more important contribution to the inclusive cross-section B — > Xs£~£+and hopefully will 
be amongst the first of the processes for which data will become available. Analysis of this process in the SM and 
in the minimal extension of the standard model have been done by many authors. Lepton polarization asymmetry 
in _B ^ K*£~ £~^wa,s first discussed by Ceng and Kao Q. In their later work they also studied SUSY effects in this 
particular decay mode 8'], which as we have already mentioned is important because it is the highest SM branching 
ratio in all the semi-leptonic decay modes. In particular Aliev et a/..lj have given the complete helicity structure of 
the amplitudes and have focused on asymmetries related to the polarization of the K* meson. Our study is more in 
the context of the simultaneous lepton polarization asymmetries and their sensitivities to various input parameters of 
the MSSM (minimal supersymmetric standard model). 

The paper is organized as follows. In the Sec. ^Iwe will present the effective Hamiltonian for the process we are 
considering, and we will write down the matrix element in terms of form factors of the B — > K* transition and then 
will give results of the partial decay rate for B K*£~£^. In Sec. IIIII we will give the analytical results of various 
polarization asymmetries. The last Sec. |^is devoted to the numerical analysis, discussion and conclusions. 



II. 



EFFECTIVE HAMILTONIAN 



The process in which we are interested {B K*£ £+) originates from the quark level transition b ^ s£ £~^. By 
integrating out the heavy degrees of freedom from the theory (MSSM here), we get the effective Hamiltonian of the 
quark level transition b 5^"^+ [11110 El 00: 



K 



4Gf 
V2 



■VtbV* 



10 



10 



(2.1) 



where Oi are current-current {i — 1,2), penguin (i = 3,..., 6), magnetic penguin (i = 7,8) and semi-leptonic 
{i = 9, 10) operators, and Ci{ii) are the corresponding Wilson coefficients renormalized at scale /i. They have been 
given in |2ll|2J|. The additional operators Qi (i = 1, . . . 10), and their Wilson coefficients are due to NHB exchange 
diagrams and are given in 0,0. 
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Neglecting the mass of the s-quark, the above effective Hamihonian gives us the foUowing matrix element: 

V Stt i_ q 

+Cq, (sPRb) {U) + Cq, (sPRb) {hd) } 



(2.2) 



where q is the momentum transfer to the lepton pair and is given as g = p_ + p+ , where p_ and p+ are the momenta 
of i~ and respectively. VtbV^g are the Cabibbo-Kobayashi-Maskawa (CKM) factors and Pl,r = (1 =F 15) I"^- In our 
analysis we will assume that we can factorize B — > K*l~ l^decay into pure leptonic and hadronic parts. ^ 

C^^^ has a perturbative part and a part which comes from the long-distance effects due to conversion of the real 



cc into the lepton pah t^i^ |g.ll8Ll2 



■res 
9 



(2.3) 



where 



^per 



Cg + -(3^3 + C4 + 3C5 + Ce) + g{m,, s)[3C, + C2 + SCs + Q + 3C5 + Cg] 
y 

-^g{l, s)[4C3 + 4C4 + 3C5 + Ce] - ^g{Q, s)[Cs + 3C4]. 



(2.4) 



The functions g{mi, s) arise from the one loop contributions of the four quark operators Oi, . . . , Og a-nd have the form 

8 4 2, 



i±V1Em. 

7(m„s) = ^iHm,) + ^ + ly^- ^(2 + y,)^\l ~ y,\ x I ''"Vi-V^^ 
y z (' y y 2 arctan 



4m? > s 



(2.5) 



where yi — Aifif/s. The non-perturbative contribution to Cg^'^ is associated with the real cc resonances in the 
intermediate states and can be parameterized by using a Breit-Wigner shape, as given in [R UslE^ : 



total 



(2.6) 



The phenomenological parameter k in the above will be taken to be 2.3 so as to reproduce the correct branching ratio 
of Br{B J/i/jK* -> K*£i) = Br{B J /il;K*)Br{J /ij} -> U). 

Using the definition of the form factors given in Eqs. lHlJ,(IX2l and we can get the amplitude governing the 
decay B K*t-l+ds^ 



M 



B^K* 



aGp 
2V2n 



VtbVZ 



{ e^.a/je*"^?"?^^ - KE + iie*.q){pK)f,F} (^7^75^) 
^iG{e*.q){U)~iH{e*.q){ej5e) 



(2.7) 



where the coefficients are 



A 



B 



S 

2mh 



2A2{s) 



+ rriK' 



c; 



(1 - rfiK') C^"T2{s) + Ai(s)(l + niK') C| 



eff 



^ There have been attempts in the Uterature to go beyond a "naive" factorization |25| . 
^ In writing this we have used q^ii'l'^i) = and 5^(^7^75^) = 2m((i'y^£). 
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1 - TTlK 



1 + rfiK' 
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1 + mK* 
E = Ai{s){l + mK-.)Cw 
2A2{s) 



1 + riiK 
2mK' 



-C 



10 



mh 
rhb 



\l + mK* s 



(2.8) 



where s — s/m\, mK* = ttik* /fnB and rhi = mi/ tub ■ From the above expression of the matrix element given in Eq. 
(|2.7() we can get the expression of the dilepton invariant mass spectra as 



dr{B ^ K*i+l-) _ Gla^m% 
ds ~ 2i07r5 



Ami 



A 



(2.9) 



where 



A = ^A(,s + 4m?)|Ap + ^^i±^(A + 12ml,:s)\B\' + l^I±^X^\C\' 



3 TO^.S 



6 rnj^,s 



-A ^^ "'f '\ s + 2m^)Re{B*C) + -{S - Arn\)\\Df 
3 mj^.s 3 

2 [A(^ + 2^,^)+^12^|,.^(g-4.h,^)] ^^^, + i-A_[A(, + 2^f) + 24^^.^,^^]|F|^ 



4Ail^#— (I + 2m,^)Re(i?*^^) + Ii:ii!M)A|GP + -A 



+2^^A[2Re(i;*i7) - (1 - mj^. - s)Re{F*H)] 



(2.10) 



and where A = A(l, s, rn^j^,) = 1 + + m^* — 2s — 2to^. — 2s Tfiif* . 



III. LEPTON POLARIZATION ASYMMETRIES 



Now we compute the lepton polarization asymmetries of both the leptons defined in the effective four fermion 
interaction of Eq. (|2.2|l . For this we define the orthogonal vectors S in the rest frame of £~ and W in the rest frame 
of for the polarization of the leptons. L, N and T correspond to the lepton being polarized along the longitudinal, 
normal and transverse directions respectively Q, IE |H 1^ - 





= (0,eL) 




= (0,ejv) 




= (0,eT) 




= (0,wi) 




EE (0,WAr) 




= (0,Wt) 



0,^ 



0, 



IP 

Pk- X p 



PK- X p_| 

(0,eArxeL) (3.1) 


' IP+I 

Q PK' X 

' Ipk- X p^ 

(0,wjvxwl) (3.2) 
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where p+, p_ and pk* are three momenta of £~ and K* respectively in the cm. frame of system. On 

boosting the above vectors defined by Eqs. I|3.2|) to the cm. frame of system, only the longitudinal vector 

will be boosted while the other two will remain unchanged. The longitudinal vectors after the boost will become 



SI ^ 



|P-I giP 
mi ' r7i<;|p_| 

|p_| EiY^ 



(3.3) 

me, me\p-\/ 

The polarization asymmetries can now be calculated using the spin projector -1(1 + 75^^) for £~ and the spin projector 
i(l + 75^) for £+. 

Equipped with the above we can now define various single lepton and double lepton polarization asymmetries. The 
single lepton polarization asymmetries are defined as Q, Q, |B>IH>I^ 



\ ds ^ ds ) \ ds 

( dTjS^, 
\^ ds 



dr{-s^,-w^) 

ds 



w^) dr{s^-w^) 

ds 



ds ds 



( dV(~S^,\ 
\ ds 

( dr(5^,-v 

\^ ds 



dr(-s^-w^ 

ds 

dr(-S:,-w, 

ds 



/ dT{S^,W^) , dT(S^-W^) 
\ ds ds 



/dr(-swM4 

\^ ds 



dr(-s^,-w^ 

ds 



(3.4) 



where the subindex x is L,N or T. denotes the polarization asymmetry of the charged lepton Along the same 
lines we can also define the double spin polarization asymmetries as Q 



( dT{S:,,Wy 

\ ds 



I _ dri-s^,Wy) 

ds 



)-( 



dT{S:,-Wy) _ dr{-s^-Wy) 

ds ds 



( rfr(s^ 

\^ dt 



Wy) dT{-S^,W, 
ds 



dT{S.,.-Wy) dT{-S:,-Wy) 
ds ds 



(3.5) 



where the subindex x, y are L^N or T . 

The expressions of the double polarization asymmetries are 



V 



LL 



Vln 



Vnl 
"Pnn 



2 (s - 4m2) 



-3^ 



.2 ^ (1 - - s)Re(B*C) + -\(\mi - s)\D, , - 
mj^tS 3 3 



2 2 [X{10mj -s) + I2srh]^, {Amj - §)] 



m,j^,s 



1 A 



6 mj^,s 



A 



[X(10mj - s) + 24:Smj^,mj]\F\- 



A 



-^^(1 - mi, - s)(10m^ - s)Re(E*F) 
3 m,j^,s 



A 



+ -r3-(s - 4m^)|Gr + -r2-s|^r + ^mi-^[2Re{E* H) - (1 - 



1 TTrhi X 
2m\, V I 



TO 



3)Re(F*i?)] 



if* 



/A 



A <^ Im(C*£:) 



(J — 4to^) 
2m£ 



Iuy{F*G) 



2rhe 



lui{C*H) 



(1 - to|-. - J) {lm{B*E) 



(1 - mj^, - s)Im{B*F) + -Im(C*F) + -^lm{B*H) - if-^I^Im{E* G) 

2 rhi rhi 



/A 



nii TT ^A(s — Amlj) 



(1 — ?fi 



2to£ 



{Re{E*H) - Re{B*G)) 



(A + 2TO^.s)Re(£;*i^) - 2sto^. {Re(B*i:') + Re{A*E)} + —X {Re(C*G) - Re{F*H)} 



-Vln 
2 X 



3 mj^. 



ml,{s-4mj)i\A\' 



\D\ 



2 m 



-2 



24 



4to^ 



mj^tTnj 



/A 



E[ 



(3.6) 



(3.7) 



(3.8) 
(3.9) 
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2m^ 



\E? - + (1 - '^i- - (l + ^ ) {Re(B*C) - Re(^*F)} 



_ 4m2)|G'p + ^s\H\^ + 6m^ {2Re(£;*iy) - (1 - m|.. - s)Re(F*ff)} 



/A 



(3.10) 



Vnt 



2X1 4m2 
3 r7i|.. 



2smlAm{A*D) - (1 - m^. - s) {Im{F*B) + lin{E*C) + 3mfIm(G*F)} 



--Im(F*C) + 2Im(S*B) + 6TO£lm(G*£;) - 3slm{G*H) 



/A 



(3.11) 



TL 



TTUlf 



y^A(J - 4m^) 



(1 - TO^,. - s) \ \E\^ + + ^ (Re(B*G) + Re(S*i/)) - ARe(£;*i^) 



1 



+2ml,s{Re{B*D) + Re{A*E)} - sX { Re(G*C) - -Re{H*F) 



/A 



Vtn = 
T'tt — 



—Vnt 
2 X 



(3.12) 
(3.13) 



3 srhl. 



m\.s{{s + Arh'l)\A\^ - (s- 4m^)|Dp} - {A(s - 2m?) - 24sm?TO^.} jSp - -X{s - Arhf)\C[' 



+ (1 - m|.. - s){s - 2rnf}Re{C*B) - (10m| - s)\E\'^ - - {(10to| - s)A - 24sTO|..m|} jFp 
+ (1 - ml, - s){l^ml ~ s)Re{E*F) + h{s ~ Amj)\G\^ - ^\H\' 



-SrfiiS {2Re{E*H) - (1 - ml, ~ sh)Re{H*F)} 



/A 



(3.14) 



where A is given in Eq. H2.10|l . 

From their definitions, Eqs. (3.1)-(3.5), polarization asymmetries relating the longitudinal (L) and transverse (T) 
spin orientations are parity odd wheras the normal one (N) is parity even. Consequently of the various double 
polarization asymmetries, Eqs. (|3.6|) - (|3.14|) . only Vln and Vtn sue parity odd. However, the basic weak interaction 
Hamiltonian is not invariant under parity transformation so that from parity symmetry considerations alone, no 
conclusion can be drawn about the vanishing or otherwise of these asymmetries. 

Since we are dealing with local Lorentz invariant theories, time reversal invariance is synonymous with CP invariance. 
In the decay process K*i^£~ , neither the initial nor the final state is an eigenstate of CP so that CP invariance 

or otherwise of the theory relate amplitudes of this process with its conjugate process B K*£~^£~ . It should be 
noted that there are terms in our matrix element which involve a triple product and thus naively have the appearance 
of a T-odd interaction. This is not correct since we are dealing with an effective Hamiltonian which includes the effect 
of strong phases which gives fake CP-violation signals even when the basic Hamiltonians are all CP-conserving. 

For the charge conjugate process the corresponding amplitudes will have their CKM factor conjugated. For 6 — s- s 
type of transition like the one considered here, the CKM phase becomes an overall phase factor since we can neglect 
the very small 6 — > it couplings. Possible CP violating phases in the CKM factor thus will not show up in any decay 
rate. Other possible sources of CP violation, for example, can come from the supersymmetry breaking parameter fi 
becoming complex. The present calculation however takes all supersymmetric breaking soft terms in the Lagrangian 
to be real so that we have effective CP-invariance of our results. The implications of these for possible measurements 
of double polarization asymmetries are remarked upon at the end of Sec. IIVI 



IV. NUMERICAL ANALYSIS, RESULTS AND DISCUSSION 



We have performed the numerical analysis of all the kinematical variables which we have presented in Sec. IIIII The 
parameters which we have used in our numerical analysis are listed in Appendix^ We have quoted our averaged 
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Standard Model values of all these variables in Table 

TABLE I; Our standard model predictions of the averaged value of the observables. 



Br(B ^ K*T-T+) 


Vll 


Pln 


Vlt 


Vnl 


Vnn 


Vnt 


Vtl 


Vtn 


Vtt 


1.29 X 10"^ 


- 0.299 


- 0.09 


- 0.329 


0.09 


-0.036 


- 0.0016 


-0.037 


0.0016 


-0.011 



We have also analyzed the effects of supersymmetry on the observables. For the numerical analysis we have 
considered MSSM, this is the simplest of the SUSY models with the least number of parameters. One of the major 
parameters of MSSM is tan/3 which is the ratio of the vev (vacuum expectation value) of the two Higgs doublets 
of MSSM. We will focus on the MSSM parameter space at large tan (3. The reason for this being that in the large 
tan [3 region of MSSM parameter space the contributions of NHB exchange becomes very important for quark level 
semi-leptonic transitions h — s- especially when final state lepton is either a muori ii or tau r. This point has 

been noted in many FCNC semi-leptonic 0| and pure dileptonic transitions [ill [l3| . Actually if we consider 
MSSM then we have to extend the set of SM Wilson coefficients, for semi-leptonic transitions we have to introduce 
two new Wilsons, namely, Cg^ and Gq^ . These coefficients come from the exchange of NHBs and are proportional to 
mbrrii tan^ /3/mh, where rrii, nib and irih are lepton, b-quark and Higgs boson mass respectively. So as we can see that 
if lepton is either fi or t and the Higgs mass is suitable then the new Wilsons (Cqj and Cq^) can have fairly large 
values. The values of Cg^ and Cg^ also depend on other MSSM parameters like chargino masses and splittings, stop 
masses and splittings etc. But as is well known these masses and splittings arc constrained by the process B Xgj 
[2?! . In our numerical analysis we will take a 95% C.L. bound [2^: 

2 X 10"'' < Br(B ^ X^7) < 4.5 X 10"'' (4.1) 

which is agreement with CLEO and ALEPH results. 

We shall now discuss the models used in our numerical analysis. The MSSM is defined on the basis of four basic 
assumptions (for a review of the MSSM refer to (i) Minimal gauge group, which is SU{3)c x SU{2)l x U{1)y 

which is the SM group also, (ii) minimal particle content, (iii) R-parity conservation, (iv) minimal set of soft SUSY 
breaking terms. If we use only these conditions then the model which is constructed is called the unconstrained MSSM 
(also called the phenomenological MSSM as one can readily study the phenomenology of it). But this sort of model 
gives rise to many phenomenological problems like FCNC, unusually large CP violation, incorrect value of Z mass etc. 
But these sorts of problems can be resolved once we make some assumptions such as all SUSY breaking parameters 
are real and hence no new source of CP violation, matrices for sfermion masses and trilinear couplings are diagonal 
which prevents tree level FCNC processes, first and second generation sfermion universality which helps us in getting 
away with the K° - K° mixing problem. 

But there is another way of solving all the problems of the unconstrained MSSM model, which is to require all the 
soft SUSY breaking parameters have a universal value at some GUT (grand unified theory) scale. If we make the 
universal values of these parameters real then even the CP violation problem is solved. This is the case in case of 
constrained MSSM and minimal supergravity (mSUGRA) models. 

Aside from the universality of all the gauge coupling constants in mSUGRA models the other conditions are: 
universality of all the scalar masses, unification of all the gaugino masses and universality of all the trilinear couplings 
at the GUT scale. With all these constraints if we impose the condition of correct electroweak symmetry breaking 
then we have another parameter which is sgn(/i)'^ and tan/3 which is the ratio of the vev of the Higgs doublets. So in 
all the mSUGRA frameworks have five parameters 

m, M, A, tan/3, sgn(/^). 



^ is the SUSY Higgs mass parameter. 
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But it is interesting to study the departure of these sorts of models. By departure we mean what would happen if 
we relax some of the above mentioned conditions of mSUGRA model. With this sort of relaxing of conditions we 
effectively introduce additional parameters in the model. One can study such relaxed models also and have reasonable 
predictions of such SUSY models if the number of new parameters introduced is not large."* There can be many options 
available; such as relaxation of universality of gaugino masses at GUT, relaxation of universality of scalar masses at 
GUT etc. 

In our analysis we will choose to relax the condition of univarsality of the scalar masses at GUT. We will assume 
non-universality of sfermionic and Higgs masses, i.e. the sfermions and Higgs have different universal masses at GUT 
scale. This sort of model we will call the rSUGRA model. With this sort of relaxation we have to introduce another 
parameter, this parameter we will take to be the mass of pseudo-scalar Higgs boson mass ruA- 

We shall now discuss the constraints put on the parameters of our models. We will consider only that region of 
parameter space which satisfies the B — > Xgj constraints given in Eq. (|4.1|l . Within the SM this decay is mediated by 
loops containing the charge 2/3 quarks and W bosons. For the set of parameters given in Appendix IbI our SM value 
of Br{B ^s7) turns out to be 3.4 x 10"'*. In SUSY theories there are additional contributions to 6 — > 57 which 
come from the chargino-stop loop, top quark and charged Higgs loop and loops involving gluino and neutralinos.^ 
Also this branching ratio constrains only the magnitude of Cj^^ . For sgn(/i) > the chargino-stop contribution 
interferes destructively with SM and charged Higgs contribution.^ The chargino stop contributions grows with tan (3 
and because of its destructive interference with the SM and charged Higgs contributions can give us a region of allowed 
parameter space. Recently there have been calculations about the NLO QCD corrections to the & ^ 57 decay rate in 
SUSY [i^l but for our work we will use the LO calculations as far as the SUSY corrections are concerned liH EB| • 

As has been emphasized in many works 0, the universality of scalar masses is not a constraint in SUGRA. To 
suppress large if" — mixing, the requirement is that all squarks should have universal mass at GUT scale. So 
that one can relax the condition of universality of scalar masses at GUT scale. This sort of model we have called 
rSUGRA. The advantage of this model arises as here we can have some handle on the Higgs boson mass and as has 
been emphasized earlier in many works the new Wilson coefficients Cq-^ and Cq^ are very sensitive to Higgs masses. 
So in this sort of model one can more easily see the dependence of various observables on the new Wilson coefficients. 

We also present the results of the average polarization asymmetries. The averaging is defined as 

_ J(3.646+0.02)Vm^^ aS 
./(3. 646+0. 02)2/m^ ds 

Although we have given the expected values of all the double polarization asymmetries with the SM in Table but 
in the graphs we have shown only those polarization asymmetries whose integrated values exceeds 0.1 either in the 
SM or in the various SUGRA models we have considered. 

In Fig. ^we have plotted the variation of differential decay rate with the scaled invariant mass of the dileptons. In 
Figs. I2l7l we have plotted the various double polarization asymmetries. In Fig. |Slwe have shown the variation of the 
branching ratio oi B ^ K*T~T'^a,s a functions of the pseudo-scalar Higgs mass in the rSUGRA model. In Fig. ^we 
have shown the variation of branching ratio as a function of tan p in the mSUGRA model. Similarly in Fig.s I10I15I 
we have shown the variation of the various integrated double polarization asymmetries as a function of the mass of 
the pseudo-scalar Higgs boson mass niA in the rSUGRA model for various values of tan/?. In Figures [161211 we have 
shown the variation of various integrated double polarization asymmetries as a function of tan (3 in the mSUGRA 
model for various values of m (the unified mass of sleptons and squarks at GUT scale) . 



Effectively this sort of model lies somewhere in between the unconstrained MSSM and th e mSUG RA model. 
^ The contribution due to the loops involving gluino and neutralinos are small as shown in I27tl30l . 

^ In our sign convention for fi it appears in the chargino mass matrix with a positive sign. 
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It is clear from the figures that several of these polarization asymmetries are sizable and that they are sensitive 
to the inclusion of the supersymmetric contributions both in regards to the magnitude and sometimes with regard 
to the sign also. The SM predictions are quite definitive; the only parameter not yet totally fixed is the mass m^, 
however, varying this within the acceptable limits does not change the values of the various asymmetries appreciably. 
Experimental observations of these polarization asymmetries will provide useful confirmatory verification of the validity 
of MSSM in rare decays of the B meson together with other experimental signatures such as single lepton polarization, 
forward-backward asymmetry etc. 

In presenting our results we have omitted showing the values of the polarization asymmetry parameters Vnt and 
Vtn, since their values are less than 0.01 and thus would be nearly impossible for observation with or without SUSY 
contributions. However, if future experiments arise with values for these which are much larger than that, it will be a 
clear indication of physics not only beyond the SM but also beyond the MSSM within the range of parameters allowed 
by other experimental constraints. 

Finally, our results pertain to the decay B K* {pi:)£'^ {p^)£^ {p_) . As discussed in the last section, for the charge 
conjugate process with the momenta unchanged, i.e. B K* (j>k)i~ {p+)i^ {p~) the polarization asymmetries {Vij) 
will be given by ±Vji, with the negative sign for Vln and Vnt and the positive sign for the others. Observations 
of these asymmetries for B and B decays would obviously need tagging of the B mesons. Observations without 
tagging with an equal number of B and B mesons would clearly produce a null value for Vln and Vnt but would 
yield value of Vll, Vnn, Vtt and (T'lt + Vtl)- The situation will change in the CKM-suppressed related process 
B pl'^l~ where because of the presence of two terms in the effective Hamiltonian with different CKM factors, 
the CKM phase would show up in the interference term and would change sign in going from this process to its 
conjugate one. Observations of asymmetries in such a process with mixtures of B and _B, as and when they become 
experimentally accessible, would provide another way of studying the CP violation through CKM phases. 
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APPENDIX A: FORM FACTORS 



The exclusive decay B K*£ ^+can be described in terms of matrix elements of the quark operators in Eq. 12. 2|) 
over meson states, which can be parameterized in terms of form factors. For B — + K*£^£^the matrix elements in 
terms of form factors of the B ^ K* transition are 00 



{K*{pk)\{V ~ A)f,\B{pB)) = -K{mB + mK')Ai{s) + i{pB + PK)t.{£* -Pb)- 



ms + nriK* 



and 



+zg^(6*.pB)^^(A3(s) - A^{s)) + e^.a.pe*y%p^j, ^^^'^ (Al) 

S TUB + fTlK" 



(i^*(pK)|s(T^.g'^(l +75)&|B(Pb)> = *e^.o/3e>gp^2Ti(s) +T2(s){e;(TO| -TO^) - [e* .pb){pb +PK)^?s 



T^{s){t\pB)\q^,-— —{pb+Pk)A (A2) 

where in the above equations pK and e^, are the four momentum and polarization vector of the K* meson respectively. 
By using the equations of motion we can get a relationship between the form factors as 

ms + niK' . , . mB ~ mK- ... ... 
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TABLE II: Form Factors for B ^ K* transition 





F(0) 


Cl 


C2 


Ai(s) 


0.337 


0.602 


0.258 


A2{S) 


0.282 


1.172 


0.567 


Ao{s) 


0.471 


1.505 


0.710 


V{s) 


0.457 


1.482 


1.015 


Ti{s) 


0.379 


1.519 


1.030 


T2{s) 


0.399 


0.517 


0.426 




0.260 


1.129 


1.128 



To get the matrix element of the scalar and pseudo-scalar currents are arrived at by multiplying Eq. (|A1|I by on 
both the sides: 

{K*{pkM^ ± l5)b\B{pB)) = -2i^{e*.q)Aois). (A4) 

mb 

For the form factors we use the results given in [s^ where we parameterize the form factors as 

F{s) = F(0) exp(cis + C2S^). (A5) 
The related parameters (ci and C2) are given in Table UTI 

APPENDIX B: INPUT PARAMETERS 

ruB = 5.26 GeV, mt = 4.8 GeV, = 1.4 GeV, 
TO^ = 0.106 GeV, = 1.77 GeV, 
myj = 80.4 GeV, = 91.19 GeV, 
VtbVtl = 0.0385, a = mK- = 0.892 GeV, 
Tb = 4.22 X 10-13 GeV, 
Gf = 1.17 X 10-5 GeV-2. 



[1] T. M. Aliev, M. K. Cakmak and M. Savci, Nucl. Phys. B 607, 305 (2001) |arXiv:hep-ph/0009133| ; T. M. Aliev and 

M. Savci, Phys. Lett. B 481, 275 (2000) |arXiv:liep-ph/0003188| . 
[2] S. Rai Choudhury, A. Gupta and N. Gaur, Phys. Rev. D 60, 115004 (1999) arXiv:hep-ph/9902355' ; S. Fukae, C. S. Kim 

and T. Yoshikawa, Phys. Rev. D 61, 074015 (2000) arXiv:he p-ph/9908229 ; T. M. Ahev, M. K. Cakmak, A. Ozpineci 

and M. Savci, Phys. Rev. D 64, 055007 (2001) |arXiv:hep-ph7 0103039 ; D. Guetta and E. Nardi, Phys. Rev. D 58, 012001 

(1998) arXiv:hep-ph/9707371 . 

[3] T. M. Ahev, M. K. Cakmak, A. Ozpineci and M. Savci, Phys. Rev. D 64, 055007 (2001) |arXiv: hep-ph /0103039| ; 

T. M. Ahev, M. Savci, A. Ozpineci and H. Koru, J. Phys. G 24, 49 (1998) arXiv:hep-ph/9705222). 
[4] S. R. Choudhury and N. Gaur, Phys. Rev. D 66, 094015 (2002) |arXiv:hep-p'h702 06128 ; G. Erkol and G. Turan, JHEP 

0202, 015 (2002) arXiv:hep-ph/0201055'. 
[5] S. Rai Choudhury N. Gaur and N. Mahajan, Phys. Rev. D 66, 054003 (2002) |arXiv:h ep-ph/0 203041| ; S. R. Choudhury 

and N. Gaur, arXiv: hep^h/0205076| E. O. Htan and G. Turan, Phys. Rev. D 61, 034010 (2000), |arXiv:hep-ph/9906502| ; 

G. Erkol and G. Turan Acta. Phys. Pol B 33, 1285, (2002) 'arXiv:hep-ph/0112115 ; G. Erkol and G. Turan, Phys. 

Rev. D 65, 094029 (2002), arXiv:hep-ph/0110017 ; T. M. Aliev, A. Ozpineci, M. Savci, Phys. Lett. B 520, 69 (2001), 

|arXiv:hep-ph/010527'9l . 

[6] F. Kriiger and L. M. Sehgal, Phys. Lett. B 380, 199 (1996), |arXiv:hep-ph/9603237| ; J. L. Hewett, Phys. Rev. D 53, 4964 
(1996), .arXiv:hep-ph/9506289| . 



11 



[7] C. Q. Ceng and C. P. Kao, Phys. Rev. D 57, 4479 (1998). 

[8] C. H. Chen and C. Q. Ceng, Phys. Rev. D 66, 094018 (2002) arXiv:hep-ph/0209352| ; C. H. Chen and C. Q. Ceng, Phys. 

Rev. D 66, 014007 (2002) arXiv:hep-ph/0205306 . 
[9] W. Bensalem, D. London, N. Sinha and R. Sinha, Phys. Rev. D 67, 034007 (2003) |arXiv:hep-ph/0209228| . 
[10] Naveen Gaur, ■arXiv:hep-ph/0305242 

[11] S. R. Choudhury and N. Gaur, Phys. Lett. B 451, 86 (1999), |arXiv:hep-ph/9810307|; A. J. Buras, P. H. Chankowski, 
J. Rosiek and L. Slawianowska, Phys. Lett. B 546, 96 ( 2002) |arXiv!hep^^pl7o207^^ J. K. Mizukoshi, X. Tata and 
Y. Wang, Phys. Rev. D 66, 115003 (2002) arXiv:hep-ph/0208078 ; T. Ibrahim and P. Nath, Phys. Rev. D 67, 016005 
(2003) arXiv:hep-ph/0208142 ; C. S. Huang and W. Liao, Phys. Lett. B 538, 301 (2002) arXiv:hep-ph/0201121, ; S. Baek, 
P. Kg and W. Y. Song, Phys. Rev. Lett. 89, 271801 (2002) arXiv:hep-ph/0205259l|. 

[12] Z. Xiong and J. M. Yang, NucL Phys. B 628, 193 (2002) arXiv:hep-ph/0105260 ; C. Bobeth, A. J. Buras, F. Kruger and 
J. Urban, NucL Phys. B 630, 87 (2002) ^arXiv:hep-ph /Ol 12305, ; C. Huang, W. Liao and Q. Yan, Phys. Rev. D 59, 011701 
( 1999) , arXiv:hep-ph/9803460 . 

[13] C. Bobeth, T. Ewerth, F. Kruger and J. Urban, Phys. Rev. D 64, 074014 (2001) arXiv:hep-ph/0104284 . 

[14] W. Skiba and J. Kalinowski, NucL Phys. B 404, 3 (1993); H. E. Logan and U. Nierste, NucL Phys. B 586, 39 (2000), 
rarXiv:hep-ph/0004139 ; Y. Dai, C. Huang and H. Huang, Phys. Lett. B 390, 257 (1997), arXiv:hep-ph/9607389 . 

[15] S. Baek, P. Ko and W. Y. Song, |arXiv:hep-ph/0208TT2l 

[16] T. M. Aliev, A. Ozpineci, M. Savci and C. Yuce, Phys. Rev. D 66, 115006 (2002) arXiv:hep-ph/0208128 ; T. M. Ahev, 
A. Ozpineci and M. Savci, Phys. Lett. B 511, 49 (2001) arXiv :hep-ph/0103261 ; T. M. Aliev, A. Ozpineci and M. Savci, 
Nucl. Phys. B 585, 275 (2000) arXiv:hep-ph/ 0002061| ; T. M. Aliev, D. A. Demir and M. Savci, Phys. Rev. D 62, 
074016 (2000) arXiv:hep-ph/9912525 ; D. A. Demir, K. A. Oliev and M. B. Voloshin Phys. Rev. D 66, 034015 (2002) 
[arXive hep-ph/0204119 ; G. Erkol and G. Turan, Nucl. Phys. B 635, 286 (2002) arXiv:hep-ph/ 0204219, . 

[17] T. M. Aliev, C. S. Kim and Y. G. Kim, Phys. Rev. D 62, 014026 (2000) |arXiv:hep-ph/991050ll 

[18] Q. S. Yan, C. S. Huang, W. Liao and S. H. Zhu, Phys. Rev. D 62, 094023 (2000) arXiv:hep-ph/0004262 . 

[19] T. M. Aliev and E. O. Iltan, Phys. Lett. B 451, 175 (1999) arXiv:hep-ph/9804458 ; T. M. Aliev, A. Ozpineci and M. Savci, 
Phys. Rev. D 56, 4260 (1997) arXiv:hep-ph/9612480 . 

[20] E. O. ntan , Int. J. Mod. Phys^A 14, 4365 (1999), [arXive |iiep^ph/9807 256 ; T. M. Aliev and M. Savci , Phys. Rev. D 
60, 014005 (1999), [arXive "hep-ph /98 1 2272 . 

[21] P. L. Cho, M. Misiak and D. Wyler, Phys. Rev. D 54, 3329 (1996) [arXi^:hep-ph/9601360. ; J. L. Hewett and J. D. Wells, 
Phys. Rev. D 55, 5549 (1997) |arXiv:hep-ph/9610323| . 

[22] D. S. Du and M. Z. Yang, Phys. Rev. D 54, 882 (1996) arXiv:hep -ph/9510267| ; T. M. Aliev, D. A. Demir, E. Iltan and 
N. K. Pak, Phys. Rev. D 54, 851 (1996) arXiv:hep-ph/9511352^. 

[23] Y. Grossman, Z. Ligeti and E. Nardi, Phys. Rev. D 55, 2768 (1997) |arXiv: hep-ph/9607473l , 

[24] B. Grinstein, M. J. Savage and M. B. Wise, Nucl. Phys. B 319, 271 (1989) ; A. J. Buras and M. Miinz, Phys. Rev. D 52, 

186 (1995) arXiv:hep-ph/9501281 . 
[25] M. Beneke, T. Feldmann, D. Seidel, Nucl. Phys. B 612, 25 (2001) | arXiv:hep-ph/0106067| ; Thorsten Feldmann, Joaquim 

Matias, JHEP 0301,074 (2003) [arXiv: hep-ph/0212158 
[26] A. Ah, T. Manuel and T. Morozumi, Phys. Lett. B 273, 505 (1991); C. S. Lim, T. Morozumi and A. I. Sanda, Phys. Lett. 

B 218, 343 (1989); N. G. Deshpande, J. Trampetic and K. Panose, Phys. Rev. D 39, 1461 (1989); P. J. O'Donnell and 

H. K. Tung, Phys. Rev. D 43, 2067 (1991) . 
[27] J. L. Lopez, D. V. Nanopoulos, X. Wang and A. Zichichi, Phys. Rev. D 51, 147 (1995) arXiv:hep-ph/9406427" ; R. Barbieri 

and G. F. Giudice, Phys. Lett. B 309, 86 (1993) arXiv:hep-ph/9303270 ; T. Goto and Y. Okada, Prog. Theor. Phys. 94, 

407 (1995) |arXiv:hep-ph/9412225| . R. Garisto and J. N. Ng, Phys. Lett. B 315, 372 (1993) arXiv:hep-ph /930730T1 . 
[28] B Physics at Tevatron : Run II & Beyond , K. Anikeev et aLarXiv : hep-ph/0201071 ; CLEO collaboration, T.E.Coan et 

aLPhys. Rev. Lett, bf 84, 5283 (2000) arXiv:hep-ex/9908022| ; ALEPH Collaboration, R.Barate et al.Phys. Lett. B 429, 

169 (1998) . 

[29] H. P. Nilles, Phys. Rept. 110, 1 (1984) ; H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 
[30] S. Bertolini, F. Borzumati, A. Masiero and G. Ridolfi, Nucl. Phys. B 353, 591 (1991). 

[31] M. Ciuchini, G. Degrassi, P. Gambino and G. F. Giudice, Nucl. Phys. B 534, 3 (1998) ^arXiv:hep-ph /9806308, . JHEP 

0012, 009 (2000) arXiv;hep-ph/0009337 . 
[32] T. Goto, Y. Okada, Y. Shimizu and M. Tanaka, Phys. Rev. D 55, 4273 (1997) arXiv:hep-ph/9609512 ; T. Goto, Y. Okada 

and Y. Shimizu, Phys. Rev. D 58, 094006 (1998) ,arXiv:hep-ph/9804294, ; J. R. Ellis, K. A. Olive and Y. Santoso, arXiv 

: |hep-ph/0204192| 



12 

[33] A. Ali, P. Ball, L. T. Handoko and G. Hiller, Phys. Rev. D 61, 074024 (2000) |arXiv:hep-ph/9910221| . 



13 




FIG. 1: Branching ratio of B ^ _ft'*r~r"'" variation with scaled invariant mass of dileptons. Parameters of mSUGRA are 
m — 200 GeV, M = 600 GeV, A — 0, tan f3 = 45 and sgn(/i) being positive. The additional parameter in rSUGRA model (the 
mass of pseudo-scalar Higgs boson) is taken to be niA = 270 GeV. 




FIG. 2: Vll variation with scaled invariant mass of dileptons. Parameters of mSUGRA are m = 200 GeV, M — 600 GeV, 
A = 0, tan/3 = 45 and sgn(/i) being positive. The additional parameter in rSUGRA model (the mass of pseudo-scalar Higgs 
boson) is taken to be rriA ~ 270 GeV. 
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FIG. 3: Vln variation with scaled invariant mass of dileptons. Parameters of mSUGRA are m = 200 GeV, M — 600 GeV, 
A — 0, tan /3 = 45 and sgn(/i) being positive. The additional parameter in rSUGRA model (the mass of pseudo-scalar Higgs 
boson) is taken to be niA = 270 GeV. 
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FIG. 4: Vlt variation with scaled invariant mass of dileptons. Parameters of mSUGRA are m = 200 GeV, M — 600 GeV, 
A — 0, tan/3 = 45 and sgn(/i) being positive. The additional parameter in rSUGRA model (the mass of pseudo-scalar Higgs 
boson) is taken to be tua = 270 GeV. 
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FIG. 5: Vnn variation with scaled invariant mass of dileptons. Parameters of mSUGRA are m = 200 GeV, M = 600 GeV, 
A — Q, tan /3 = 45 and sgn(/i) being positive. The additional parameter in rSUGRA model (the mass of pseudo-scalar Higgs 
boson) is taken to be ruA = 270 GeV. 




FIG. 6: Vtl variation with scaled invariant mass of dileptons. Parameters of mSUGRA are m = 200 GeV, M = 600 GeV, 
A = Q, tan /3 = 45 and sgn(/i) being positive. The additional parameter in rSUGRA model (the mass of pseudo-scalar Higgs 
boson) is taken to be uia = 270 GeV. 
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FIG. 7: Vtt variation with scaled invariant mass of dileptons. Parameters of mSUGRA are m = 200 GeV, M — 600 GeV, 
A — 0, tan /3 = 45 and sgn(/x) being positive. The additional parameter in rSUGRA model (the mass of pseudo-scalar Higgs 
boson) is taken to be niA = 270 GeV. 
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GeV, M = 450 GeV, A = 0. 
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FIG. 12: (Vlt) variation with m,4 (in GeV) for various values of tan f) in rSUGRA model other model parameters are m = 200 
GeV, M = 450 GeV, A = Q. 
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FIG. 14: (Vtl) variation with niA (in GeV) for various values of tan /3 in rSUGRA model other model parameters are m = 200 
GeV, M = 450 GeV, A = Q. 



20 



-0.2 



V 




-0.4 



-0.6 



200 



300 

(GeV) 



400 



FIG. 15: (Vtt) variation with (in GeV) for various values of tan /? in rSUGRA model other model parameters are m — 200 
GeV, M = 450 GeV, A = 0. 



0.4 
0.3 
0.2 
0.1 

A 

c^- 

V 

-0.1 
-0.2 
-0.3 



-0.4, 





' 1 


1 1 1 1 1 




SM 


^ — 

y 




-- m = 300 
■-■ m = 500 


y 




■-■ m = 700 


— 
















































I.I.I. 



30 



35 



40 

tanP 



45 



50 



FIG. 16: {Vll) variation with tan/3 for various sets of m in mSUGRA model. Other model parameters are M = 500 GeV, 
A = Q. 
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FIG. 17: {Vln) variation with tan/3 for various sets of m in mSUGRA model. Other model parameters are M — 500 GeV, 
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FIG. 18; (Vlt) variation with tan/3 for various sets of m in mSUGRA model. Other model parameters are M = 500 GeV, 
yl = 0. 
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FIG. 19: {Vnn) variation with tan/3 for various sets of m in mSUGRA model. Other model parameters are M — 500 GeV, 
A = 0. 
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FIG. 20; {Ptl) variation with tan/3 for various sets of m in mSUGRA model. Other model parameters are M = 500 GeV, 
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FIG. 21: {Vtt) variation with tan/3 for various sets of m in mSUGRA model. Other model parameters are M = 500 GeV, 



